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We show that conventional asymmetric chiral random ma- 
trix models (ChRMM), with a gaussian distribution in the 
asymmetry, provide for a screening of the topological charge 
and a resolution of the (7(1) problem in the unquenched ap- 
proximation. Our exact results to order 1/N are in agreement 
with numerical estimates using large ensembles of asymmetric 
ChRMM with gaussian distributions. 



I. INTRODUCTION 

In QCD the axial U(l) symmetry is not apparent in 
the spectrum, although a non-vanishing quark conden- 
sate suggests the existence of a ninth Goldstone boson. 
This is the U(l) problem. The resolution of this appar- 
ent problem is believed to follow from the chiral anomaly 
jij. 't Hooft has suggested a specific mechanism using 
instantons ||]. Witten has proposed a resolution in the 
context of the large N c (number of colors) limit ||, an 
idea that was interpreted by Veneziano in terms of vector 
ghosts ||. 

In all existing scenarios for the resolution of the U(l) 
problem, it is crucial that the topological susceptibility 
is nonzero in both the quenched and unquenched (with 
massive quarks) approximation. Lattice simulations ap- 
pear to support this assumption both for quenched g] 
and unquenched although the idea may be at odd 
with translational invariance and current identities |Q. 
This notwithstanding, it was argued by few that a sta- 
tistical ensemble of topological charges may yield to the 
screening of the bare topological susceptibility and the 
resolution of the U(l) problem ||. This scenario will be 
considered in the context of standard ChRMM |9|dT^] . 

Standard ChRMM follow from the constant mode sec- 
tor of the instanton liquid model Jl3[-^5[, regarded as a 
statistical ensemble of topological zero modes. As mod- 
els, they offer a minimal framework for discussing the 
interplay between (chiral) symmetry, typical scales and 
the thermodynamical limit. In these models gaussian 
randomness is enough to cause the spontaneous breaking 
of chiral symmetry (hermitean matrices) or the sponta- 
neous breaking of holomorphic symmetry (nonhermitean 



matrices) in the thermodynamical limit. They have also 
been used to investigate issues related to the universal- 
ity of "noise-fluctuations" in Dirac spectra both in the 
microscopic || and macroscopic ]l6| limit. 

In so far, most of the analyses carried in the context of 
standard ChRMM have been conducted using standard 
and symmetric ChRMM. In the context of the instanton 
liquid model this means that the number of instantons 
and the number of antiinstantons is fixed. In this case 
there is no resolution of the U(l) problem. In this let- 
ter, we will consider a statistical ensemble of asymmetric 
ChRMM with a gaussian distribution for the asymme- 
try. The physical motivation for this model stems from 
the coarse grained version of the instanton liquid model 
p9[ . In section 2, we will discuss basic QCD Ward iden- 
tities emphasizing the role of a non- vanishing topological 
susceptibility in the resolution of the U(l) problem. In 
section 3, we streamline the phenomenological arguments 
for standard but asymmetric ChRMM with a gaussian 
distribution for the matrix asymmetry. In section 4, we 
use a diagrammatic analysis based on a 1/7V expansion 
to solve the U(l) problem. Our analytical results for the 
unquenched topological, scalar and pseudoscalar suscep- 
tibilities are in agreement with numerical estimates using 
large asymmetric matrices. Our conclusions are in sec- 
tion 5. We elaborate on some derivations in the Appen- 
dices. Throughout, we will use four-dimensional argu- 
ments in Minkowski space. The transcription to ChRMM 
will be done through their Euclidean four-dimensional 
analogue. 



II. QCD WARD IDENTITY 

The topological term in the QCD Lagrangian 9 3 is a 
total divergence 



32^2 



(1) 



where is the Loos-Chern-Simons current. For Nf 
flavors of quarks with current masses m/, the gauge- 
invariant flavor axial-singlet current j^" v is anomalous, 



d»j™ = 2N f E + 2j2m f ^ f i^f 
f 



(2) 
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In a 9 state, the expectation value of (g) implies 

N, 

U = 2N f {9\E\9) +2j2(#\™ f ipfi'rsi>f\6) (3) 
/ 

by translational invariance. In particular, the topological 
susceptibility is given by 



= ^W 1 = -^7 E ^\m^ f i^f\6) (4) 



The absence of a physical massless £7(1) boson for 
quarks of equal masses m/ = m, gives [pj 



= J d^(9\T*f^(x)^ l5 m\0) 

_ A? 

= -2i<0|^|0) - 2i^x(0) 
m 

+2m / d 4 x(6»|T* ; 0i7 5 V'(x) ^75^(0)161) (5) 



where we have made use of (||,[|) . T* is the covariantized 
T-product. For 9 = 0, it follows from (§) that 



the (7(1) anomaly, have been put forward by several au- 
thors (20|. In these effective models it is important that 
point-like pseudoscalars are coupled to point-like glue- 
balls to achieve consistency with the QCD anomalies (ax- 
ial and scale anomaly). 

Such effective models arise naturally from microscopic 
descriptions of the QCD vacuum using a random en- 
semble of instantons and antiinstanto ns p 3H15i. In the 
"coarse grained instanton model" J15|fL9[ 
"glueball" fields are identified with 



the effective 



32tt 2 
32tt 



G-G(x)-> (n+ + n-)(s) 
-G-G{x)^ (n + -n-)(x) 



(8) 
(9) 



where n±(x) is the density of instantons (+) and anti- 
instantons (— ). The U(l) anomaly is saturated by as- 
suming that the distribution in the number difference 
(ri + — n_) (susceptibility of the quenched vacuum) is 
gaussian. The scale anomaly is also saturated by assum- 
ing that the distribution in the number density (n++n_) 
(compressibility of the quenched vacuum) follows from a 
logarithmic ensemble Jl9|,p0[. 



N] 



Iffl t , , __ 

iXtop = + TFI I d 4 x(T*ipij 5 ip(x) ipi-f 5 tp{0)) (6) 



where we have set Xtop = x(0), in agreement with earlier 
results |L7|,[l8|. The resolution of the U(l) problem stems 
from the observation that for small m, the absence of a 
U(l) Goldstone mode requires that \to P — —m(ipilj)/Nj 
to order C(m 2 ). 

A rerun of the above argument for the SU(Nf) cur- 
rents and densities yields the flavor non-singlet relation 



A. Symmetric ChRMM 

Standard ChRMM are schematic versions of current 
descriptions of the instanton liquid model Jl3|-p^|. One 
essentially truncates the QCD partition function to the 
space of instanton-antiinstanton zero modes and further- 
more assumes that the fermionic overlaps are randomly 
distributed with a Gaussian measure. Specifically 



Z 



det N, Q 



(10) 



= ~<V [A/,A, 7 ] + ^> 



where 



+ra / d i x(T*Tjjij 5 \itp(x) ipi^Xjipfi)} 



(7) 



which shows that for small m there should be a multiplet 
of Nf — 1 Goldstone modes. Contrasting (0) with (|7]) 
shows the importance of a non- vanishing topological sus- 
ceptibility in the resolution of the U{1) problem. It was 
originally argued by Witten || that further consistency 
with large N c arguments requires that the quenched topo- 
logical susceptibility x* be nonzero as well. Current lat- 
tice simulations seem to support this conjecture |J, al- 
though there may be subtleties as we indicated above 



III. CHRMM 

A number of effective models aimed at describing the 
long wavelength physics of the QCD vacuum including 



Q 



inif A* 
A irtif 



(11) 



and A is an N X N matrix distributed with a weight 



P{A) 



-JVEtrAA* 



(12) 



where E is a fixed scale related to the quark conden- 
sate through < ip'ip >= l/VE. For simplicity E = 1 in 
the arguments to follow. The large N limit is directly 
linked to the thermodynamic limit through the choice 
N/Vi = 1. In other words the number of zero modes is 
made commensurate to the four volume V4, to insure a 
non-vanishing chiral condensate in the thermodynamical 
limit. 

One can easily check that for such ensembles, the topo- 
logical susceptibility is zero. This is due to the fact that 
in ( |io| ) we have used symmetric matrices, i.e. in the 
instanton liquid model, we have set the number of in- 
stantons equal to the number of antiinstantons. 
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B. Asymmetric ChRMM 



A. The Rules and Resolvent 



The above constraint can be relaxed following 
HIPP! through 



+ ,n- 



(13) 



• Feynman Rules 

To assess (|l7|) we use a 1/N expansion. The Feynman 
rules associated to ( |l5| ) are shown in Fig. [|. The bare 
quark propagator is l/(— im) and the bare gluon propa- 
gator is 



where \ = n + - 



is the difference between the number 



of instantons and antiinstantons and Z n+tTl _ is as before _ (A a A°\ = 

©, but with the A's now rectangular and n+ x n_ 
complex valued. The star in ( |l3| ) indicates that the sum 
is restricted to n + + n_ = 2N with N eventually going to 
infinity. The issue related to the scale anomaly (for 
our case an unrestricted sum in dig))) will be discussed 
elsewhere. The topological susceptibility in the limit of 
infinite quark mass (quenched) is in this case just 
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Xtop 



N 



(xx) = X* 



(14) 



Although ( p"3] ) describes a statistical ensemble of point- 
like topological defects, that is zero correlation length, 
the latters are spread over a gaussian of width \/N\*, 
resulting into a finite correlation length or topological 
susceptibility in the thermodynamical limit. In the next 
section we evaluate the pseudoscalar and topological sus- 
ceptibilities for the ensemble (|lj) in the context of a 1 /N 
expansion. 



IV. U{1) PROBLEM IN CHRMM 

A measure of the axial-singlet charge in the QCD vac- 
uum in the presence of light quarks is given by the pseu- 
doscalar axial-singlet susceptibility 



(15) 



A resolution of the U(l) problem requires that x ps re- 
mains finite as V4 — > 00 followed by m — > 0. 

To help investigate ( |l5| ) in the context of standard 
ChRMM, we recast ( |l3| ) in the form of a partition func- 
tion over Grassmann variables 

Z = / dtp^dipdAe ^ N f (16) 



[-(trQ-SCTH)] 
[(tr(Q- 1 75 )tr(Q- 1 75 )) c 



(17) 



The Euclidean analogue of (|T|) in ChRMM is 
1 

Xps - 

1 

+ N 

where all averages are carried using Jl6| ) that is (13). 
A naive expectation based upon large N counting rules 
suggests that the second term in (17) is subleading in 
the thermodynamic limit. As we will show below, this 
expectation is not born out by calculations. 



Here a, b, c and d run from 1 to 2N, ® denotes ten- 
sor product and T is a flavor bearing matrix with flavor 
indices /,;. For Nf = 3 

- *?A - l^&oYfl + Mj[Az]g) (19) 

following a standard decomposition. In terms of®, the 
bare gluon propagator may be rewritten in the form 



V = — ((lAj) (8 (lAj) - (75A/) ® (75A7)) (20) 

with / running from to 8. We use the normalization 
trA/Aj = 25 1 j. The flavor decomposition stemming from 
( pp| ) is reminiscent of the flavor interaction in a schematic 
quark-meson effective theory of the Nambu- Jona-Lasinio 
type. 

_J_ 

-im A 75 



I+0+]QlT 




FIG. 1. Feynman rules. In the upper part: the bare quark 
propagator l/(— im), the dressed one (Q), the 75 vertex, the 
bare \ propagator and the bare gluon propagator T>. In the 
lower part the dressing of the gluon propagator is shown, the 
last equality representing the Schwinger-Dyson equation. 

• Resolvent 

For a single flavor we define the resolvent 



G(z) 




1 



im 

A im 



(21) 



where averaging is now over rectangular matrices. For 
our purposes we will need the resolvent calculated at z = 
only. For fixed n± , the result of the averaging over A 
yields 



S = gl + 9515 



(22) 
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with Q — Q{0) and 



^— m 2 — 2x 2 + m\J 4 + m 2 + 
m(l — a; 2 ) 



55 



i I m 2 + 2 - m-i/ 4 + m? + ^4 



m(l — x 2 ) 



and 



X 

2N 



2N 



(23) 



(24) 



Both 1 and 75 in are 2N x 2N valued. The explicit 
form of the resolvent g with m = —iz in the chiral limit, 
can be checked by other methods |2j,|2^]|^|. Both g and 
<75 are flavor additive. We observe that 



tvQ = g + xg 5 = - I -m 



^ (25) 



and 



tr7 5 C/ = xg + g 5 = 



IX 



(26) 



For x = (square matrices) the discontinuity of ( B5h 
i/+(A,a; = 0) = --Im tr£(m = -iA + 0) (27) 

7T 

is just Wigner's semicircle for the quark spectral distri- 
bution. The argument of Q in ( p7| ) has been appended 
for clarity. For a; ^ (asymmetric matrices) it is the re- 
solvent in a configuration of unequal topological charges 

u + (X,x) = \x\5{\) + ^^(A 2 - A 2 _)(A 2 + - A 2 ) (28) 



with Aj_ = 2±2Vl — x 2 . The delta function at the origin, 
reflects on the number of unpaired topological charges. 
Configurations with x 7^ do not break spontaneously 
chiral symmetry, owing to the occurrence of a gap at 
the origin. The spontaneous breaking is triggered by the 
neutral topological configurations with x = 0. 

We observe that the discontinuity in (|2^) for fixed x 

v_(X,x) = -ilm tr(7 5 <7)(m = -iA + 0) = x5(X) (29) 

is a direct measurement of the difference in the spectral 
distribution between left-handed A~ and right handed 
A+ quark zero-modes. The result ( p9| ) reflects on the 
Atiyah-Singer index theorem. 



In the large N limit, the fluctuations in \ are of or- 
der \/~N, and a statistical averaging over an ensemble of 
asymmetric matrices yields x = on the average. A 
convenient way of interpreting the averaging over \ dia- 
grammatically is to consider the partition function jl^ ) 
as describing the interaction of random matrices with a 
field x whose propagator is (xx) — ^X* ( see Fig- 0)- 





r 2 

FIG. 2. 'Trace' contribution (ri + I^) to the pseudoscalar 
susceptibility. 

• Gluon Ladder 

In the next step we consider the resummation of the 
gluon ladder. If the resolvent is diagonal]^ 

G = (-i)(a lX + a 8 lX s ) (30) 

and the resummed planar propagator is decomposed as 

1 



AN 



(DJ <g> (1A/) - Df ® ( 75 A/)) 



(31) 



in analogy with (|2^) , then the resummation of the rain- 
bow diagrams gives 



Df 



Ih + ^(tr\ I g\ J g)D s j 



(32) 
(33) 



This resummation is shown in the form of a Schwinger- 
Dyson equation in Fig |IJ The final result is 



V 



resurn m cd 



1 

2N 



1 



2 + d 



1A/ <g> lAj 



1 



2-d 



1. j 

75A/ <8) 7 5 Aj 



where d is 9 x 9 matrix valued 

du = -tiXiGXjG 



(34) 



(35) 



The denominators 2 ± d in the pseudoscalar channel can 
be interpreted within the effective quark-meson analogy 
as generalized GMOR (Gell-Mann-Oakes-Renner) rela- 
tions. Using these results we may now evaluate the pseu- 
doscalar susceptibility. 



Throughout m u = mj 7^ m s . 



1 Note that a different measure was used in |25| 
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B. The Puzzle 

Now consider the pseudoscalar axial-singlet suscepti- 
bility in ChRMM 



(■0 t A O 7 5 V''0 t Ao75V'}c 



(36) 



Let us first look at the diagrams corresponding to the first 
line ("trace term") in (17). The graphs that contribute 
to the 'trace' part are shown in Fig. [2| The "trace-trace" 
term, corresponding to second line of (0), will be con- 
sidered later. We remark that main difference between 
the singlet and nonsinglet pseudoscalar correlators is the 
second line of (|T7j). 

The first graph in Fig. || is given by 



Ti = tr£7 5 Ao<?75A = g tr £ 2 



(37) 



The prefactor of 2/3 follows from our choice of Ao = 
\/2/3. Note that we may safely set x — here. The 
second graph in Fig. || is given by 



r 2 = tv(g l5 x g)v r 



esumrne' 



(38) 



A simpler way to evaluate this graph is to use the identity 



Note that the bare "quark lines" carry only a factor 1 
(amputated propagators) . This identity comes as a result 
of contracting a bare gluonic propagator (|20| ) with a 75 Ao 
vertex. Using the definition of the resummed propaga- 
tor (Schwinger-Dyson equation) we arrive at the useful 
identity 






Applying the above identity twice, yields 

We remark that in the sum T1+T2 the Ti term cancels 
against the first one in the above line and the result for 
the sum is simply the sum of the two remaining diagrams 
in the line above, that is 



Ti + T 2 = +4A + (2N • 2) 



-1 
2N 



1 



2-d 



(27V ■ 2) 



(Ml 



4 A' 11-2 

N f 



1 



2-d 



E 



m i (y / A + mf + mi) 



(39) 



The diagrams of Fig. 2 yield a Goldstone pole in the 
pseudoscalar axial-singlet correlator, hence the a priori 
puzzle. 



Before we proceed to solve the puzzle, we note the re- 
markable cancelation in the sum Ti +r2. The self-energy 
insertions on the two-quark lines and the one-gluon ex- 
change diagram, conspire to remove the "constituent" 
quark cut, leading to a Goldstone pole. For the non- 
singlet channels this is the mechanism by which the Gold- 
stone modes emerge in a correlation function that does 
not abide by confinement. This cancelation is reminis- 
cent of the one taking place in two-dimensional QCD in 
the planar approximation p6fl . So it appears that the 
removal of the "constituent" quark cut is a simple con- 
sequence of a consistent treatment of a Ward identity 
(vertex and self-energy) without the need of confinement. 



C. The Solution 

In the previous section we have neglected the effects 
due to the asymmetry in n± as they appear to be 
subleading in 1/iV, along with the quark loops. In a 
way, we derived a quenched result. However, the oc- 
currence of strong infrared sensitive terms of the form 
x/m ~ in the diagrammatic expansion requires 

that we resum them, prior to taking N — > 00 for finite m. 
We show below, that the infrared sensitive terms upset 
naive power counting in the axial singlet channel. They 
are just a manifestation of a screening phenomenon in a 
statistical ensemble made of negative and positive topo- 
logical charges. Technically the necessity of rederiving 
1 /N counting rules stems from the fact that we have in- 
troduced here a non-large N ingredient — the x field. 

The new diagrams that can contribute to our above 
analysis are the formerly disconnected quark loops which 
now interact via \ propagators. Looking at a few dia- 
grams we conclude that 

• the x lines cannot end on the same quark loop. 
This gives a subleading contribution in 1/N. 

• when considering the Green's function ( a 1-point 
correlator ) we get no additional graphs. 



• two point correlators of the type 



(40) 



get contributions on the endpoints from graphs pro- 
portional to x interacting through a dressed \ prop- 
agator. 

• the dressed x propagator gets contributions from 
the part of a closed quark loop which is propor- 



tional to 



x 



With the above observations, we now proceed to solve 
the puzzle. 
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D. Screening 



The quark loop with 75 Aq insertion gives 



Here we isolate the contribution of the closed quark 
loop r which is of order \ 2 ■ F° r that consider a general 
Feynman graph made up of undressed quark and gluon 
propagators. A factor proportional to x can only appear 
from a summation in a quark loop (with the 75 part of 
a propagator contributing) . Therefore we must have two 
distinguished quark loops in the graph T. Now we may 
reduce the parts of T lying 'at the edges' and between the 
distinguished loops. The former transform into a trace 
of the (symmetric) Green's function, while the latter one 
becomes the resummed gluon propagator (|I|). We must 
yet ensure that each loop picks up a factor of x- For 
this we must take only the 75 channel of the resummed 
propagator. In the end we obtain: 



- (-l)tr(£ 75 A/) 



1 

2N 



1 



tr(£ 7 5Aj) (41) 



The first minus sign comes from the fermion loop, while 
the second one from the resummed gluon propagator. 
This gives (recall eq. (|30|)) 



X 2 ci 



1 

2N 



1 



2-d 



cj 



i.j 



with ci = itrQXi. Resummation gives 

1 

Xtop 



1 




1 




x» 




2-d 


Cj 

IJ 



Using the definitions for cj and du we obtain 

N f 

E 



1 



Xtop 



1 

X* 



1 



(42) 



(43) 



(44) 



For Nf = (quenched) the topological susceptibility is 
Xtop = X*> as it should. For Nf 7^ (unquenched) the 
topological susceptibility is screened by the quark loops 
and vanishes in the chiral limit |p^| . The alternative, 
non-diagrammatic proof of equality ([5]) is presented in 
Appendix A. 



E. Pseudoscalar susceptibility 

Given the above observations, we conclude that our 
previous calculation for Ti and T2 is correct. However 
our neglect of the term from the second line in ([It]) is 
not correct. Although, this term cannot receive a con- 
tribution from gluon insertions as they are suppressed in 
1/N, they may have a contribution through a fluctuation 
in the size of the quark matrices Q. This is essentially 
an exchange of a x ~ VN field as shown in Fig. |^. 



FIG. 3. 'Trace-trace' contribution (Fs) to the pseudoscalar 
susceptibility. 



N, 



tr S75A0 = 2Nd -IJ2—X + 0(x 2 ) 



3=0 1 



(45) 



so the relevant part is just \/2/3 ■ -^rX- The full graph 
gives 



(46) 




The susceptibility is then the sum of the two graphs of 
Fig. H and the graph of Fig. ||, that is 



2N 

ri+r 2 -r 3 = — 



where we used the notation 



mi(y / 4 + rof + rm) 
For equal quark masses we get 



(47) 



(48) 



2N 

Xps ^~ 



2N 



-12S 



9/m 2 

x 

125 _ 4 . 9^2 



35 



^7 + 3 * 



(49) 



Since our choice of Ao = \/2/3 and 75 is 2A^ x 2A^ valued 
with n++n- = 2N fixed, the prefactor of 2N/3 is natural 
in (^9j). The properly normalized susceptibility is 

3 3 

Xps = l^Xps = ^( r 3- r i ~ r 2) (5°) 

which is finite in the large N limit. 

In the chiral limit S ~ l/2m and the leading singular- 
ities in 1/m 2 cancel in the numerator. The pseudoscalar 
susceptibility Xps is finite in the chiral limit, thereby solv- 
ing the U(X) problem. The a priori Goldstone pole in 
Ti + T 2 was removed (screened) by the order V^V fluc- 
tuations in the size of the matrices. Although our result 
was derived with a gaussian weight in A, we believe the 
latter to be a fixed point as discussed in |27|]. 



F. Scalar Susceptibility 

Similar arguments for the scalar susceptibility 
1 



(51) 



translates in ChRMM to 



G 



Xs = ^ [-(trQ-^- 1 )] 
+ ^ [(trCQ-^trCQ- 1 ))^ 



(52) 



ignoring (ultraviolet sensitive) contributions from con- 
tinuum states. Similar correlators using also results from 
random matrix theory where also discussed in [H . 

Since the scalar susceptibility is proportional to the 
derivative of the Green's function with respect to the 
mass, we can immediately write down the result for (|52~ 
in the form 



We have checked that our numerical results are in 
agreement with the Ward identity (^J) when translated 
to ChRMM. Specifically 



Xtop - 



N 2Xps 



2mi 



N 



2 trQ 

f 



(54) 



The proof of the above identity is included in Ap- 
pendix B. Figures [I] and || (lower right) show the Ward 
identity (54). The numerical points represent the l.h.s. 
of the Ward identity, the solid line is calculated from the 
analytical prediction for the r.h.s. using (|25|) with x = 0. 



Xs = N f 1 - 



V4- 



(53) 



This result is only qualitative, since we have ignored the 
effects of the scale anomaly. Also for small values of m, 
( p3| ) does not show terms of the form lnm expected from 
the exchange of two Goldstone modes. In ChRMM these 
exchanges are 1/N suppressed. 



G. Numerical Results 

To check on the validity of the above arguments, we 
have carried out direct numerical calculations using a 
gaussian distributed ensemble of rectangular matrices. 
Figures |and | show the results of the numerical sim- 
ulations for the topological, pseudoscalar and scalar sus- 
ceptibilities, with one and three flavors respectively. The 
solid lines follow from (0), ([l9|), and (|5^) respectively. 
Our 1/N analysis is confirmed, except for very small 
values of m, where finite 1/N size effects are noted. 
Throughout x* was se t to 1. 
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FIG. 4. Normalized topological (upper left), pseudoscalar 
(upper right), scalar (lower left) susceptibilities, and Ward 
identity (lower right) for Nf = 1. The numerical simulations 
where carried for fixed 2N — n+ + n_ = 60, 120, 200 and 
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FIG. 5. Same as Fig. | but for N f = 3. 



V. CONCLUSIONS 

We have shown that an ensemble of asymmetric but 
standard ChRMM with an assumed gaussian distribution 
for the asymmetry, yields a screened topological suscepti- 
bility and solves the U(l) problem. We believe this result 
to be new in the context of ChRMM. Using diagram- 
matic techniques in the context of the 1/N approxima- 
tion, we have shown the fermion determinant contributes 
importantly in the axial-singlet channel through infrared 
sensitive terms, but otherwise is quenched in most other 
channels in the thermodynamical limit. We also believe 
that this result is new in the context of ChRMM. Some 
of these results are generic and illustrative of the mecha- 
nisms at work in more realistic models of the QCD vac- 
uum, such as the instanton liquid model. It would be 
interesting to see how the present arguments extend to 
nonstandard ChRMM |^| , and how they are modified at 
high temperature and finite N in comparison to current 
lattice simulations ffl). Also, our analysis allows for an 
assessment of the 9 vacua in ChRMM, as we will discuss 
next. 



< (n+ — n_ 
result. 
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Appendices 

A. Screened Topological Susceptibility 

In this Appendix, we will outline an alternative deriva- 
tion in favor of the screening of the topological charge in 
the present model. The connected 'vacuum' diagrams 
(prior to averaging over x) ar e given by logZ(x). For a 
fixed asymmetry x = x/2A = {n + — n_)/2A we have 



d m \ogZ( X ) = -i{ipH) = -2NitiQ 
= NN f -to + \ A + to 2 



Ax 2 



(55) 



The contribution of order x 2 follows by expanding the 
right hand side 



(d m logZ( X )h=NN f - 



2x 2 



TO 2 V '4 + TO 2 

and integrate it with respect to to to obtain 



(56) 



(logZ( X )) 2 = NN f (- Vm * +4 +co nst) x 2 (57) 
\ 2to I 

The constant is set by the requirement that for to — > oo 
we should obtain zero. Hence 



(\ogZ( X )h = NN f \- 
= —2NNf 



\Jm 2 + A 1 
2~n~i + 2 

1 



to(-\/4 + to 2 + m) 
1 



2N to(V 4 + to 2 + to) 
The full partition function is is given by 



• X 



dxe 



'Z(x)- 



Inserting (|o|) into ( j6l| ) yields 

Z = f dx e-H^ +N >- 
in agreement with (H). 



(V4+m2+m) 



(58) 
(59) 
(60) 

(61) 
(62) 



B. Ward identity in ChRMM 

The Ward identity (Q) can be derived from the par- 
tition function for finite angle and massive quarks. 
Specifically 



Z[B] = (dct 



Nj 



melNf A L )) = (det^Q 9 ) (63) 
A ime N f / 



Using the derivatives at 9 = 

in 



deQe 



N f 



75 



and 



we get 



N X 



top 



Since (ip^ip) = 2NtiG, then 
2to 



Xtop - 

in agreement with (|54|) 



^ TO _ 

Nj N 2 } Xps 



(64) 



(65) 



(66) 



(67) 
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